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We study the spontaneous emission of an excited atom close to an optical nanofiber and the 
resulting scattering forces. For a suitably chosen orientation of the atomic dipole, the spontaneous 
emission pattern becomes asymmetric and a resonant Casimir-Polder force parallel to the fiber axis 
arises. For a simple model case, we show that the such a lateral force is due to the interaction of the 
circularly oscillating atomic dipole moment with its image inside the material. With the Casimir- 
Polder energy being constant in the lateral direction, the predicted lateral force does not derive 
from a potential in the usual way. Our results have implications for optical force measurements on 
a substrate as well as for laser cooling of atoms in nanophotonic traps. 
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Electromagnetic fields close to a surface and their in¬ 
teraction with particles are of fundamental interest. Re¬ 
search in this area covers, e.g., dispersive interactions 
such as Casimir and Casimir-Polder (CP) forces [T|, as 
well as forces that arise from the scattering of external 
light fields. Recently, optical forces which act perpen¬ 
dicularly to the propagation direction of an excitation 
light field attracted increasing interest [2] - In particular, 
it has been predicted that evanescent fields exert lateral 
forces and torques on Mie particles [3] . Even stronger lat¬ 
eral forces are expected for chiral particles in evanescent 
fields g]. 

Recent experimental research on emitters close to sur¬ 
faces have revealed that suitably excited particles can be 
used to realize strongly directional excitation of guided 
modes mm- For example, when a gold nanoparticle on 
the surface of an optical nanohber scatters the light of an 
external circularly polarized laser beam, the coupling into 
counter-propagating guided modes of the nanofiber can 
exceed a ratio of 40:1 m • Such an asymmetric scattering 
is independent of the excitation process and is only gov¬ 
erned by the polarization of the emitted light. The con¬ 
servation of total momentum in the system in conjunc¬ 
tion with the asymmetric emission suggests the existence 
of a force on the scatterer that is parallel to the waveg¬ 
uide axis. However, the coupling to radiative modes has 
to be taken into account as well [T3j . The search for such 
a scattering force is within reach of current cold atom 
experiments, in which asymmetric excitation of guided 
modes has already been observed [9]. 

Lateral forces have also been discussed within the con¬ 
text of dispersion interactions where they arise even in 
the absence of external fields. Lateral CP forces are typ¬ 
ically achieved by breaking the translational invariance 
of the surface via periodic corrugations ehhh] or disor¬ 
der [IS]. Moreover, lateral Casimir forces between two 


periodically structured surfaces have been proposed 
f2Tj and measured [22], and have been put forward as a 
means to realize contactless force transmission [25] . 

Here, we propose to exploit the directional sponta¬ 
neous emission by an atom near a nanofiber to realize a 
translationally invariant lateral CP force. The envisioned 
physical situation is illustrated in Fig. [lja) and an illus¬ 
tration of an emission pattern is shown in Fig. [ijb). A 
cesium atom is located at a position ta = xax at a dis¬ 
tance cIa = xa — R from the surface of a fused-silica fiber 
of radius R = 250 nm. The atom is initially prepared 
in an excited state |1) = |6 2 P 3 / 2 , F' = 5, M' F =5) . Here, 
the y coordinate axis is chosen as the quantization axis. 
The only available decay channel is to the ground state 
|0) = |6 2 Si/ 2 , F = 4, M f =4), such that the decay of the 
atom leads to the emission of a <r + -polarized photon. The 
transition has a free-space wavelength of Aio = 852 nm, 
and corresponding wavenumber fcio = 27 t/Aio and fre¬ 
quency wio = cfcio- The transition dipole matrix element 
is dm = 1.9 x 10 -29 Cm x (ia: + z) pSi] , 

We use two methods to calculate the emission rates 
and Casimir-Polder force. The first method allows 
an intuitive interpretation and is based on solving the 
Schrodinger equation for the atom-field interaction for 
a suitable mode decomposition of the electric field op¬ 
erator [25 . The second method uses the more general 
Green’s tensor formalism where, for example, absorption 
of the nanofiber is easily included. 

The total spontaneous decay rate of an atom can be 
expressed in terms of the Green’s tensor G(r, r',u) as [T 

r {r A ) = (2 n 0 /h)ul 0 d 10 ■ ImG(r A ,r A ,wi 0 ) • doi- (1) 

A lateral force on the atom may result from an unbal¬ 
anced spontaneous emission rate into the +z and —z half 
spaces. To obtain information about the directionality of 
the spontaneous emission, we decompose the electric field 
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FIG. 1. (a) Sketch of the situation under consideration. An 
atom (yellow sphere) is located at ta = (xa, 0, 0), in close 
proximity to a nanofiber, modeled as a dielectric cylinder in 
free space. Calculations are performed for the cesium D 2 line 
optical transition at a wavelength of A 10 = 852 nm. The 
atom is initially in the excited state 11), from where it can 
only decay by emitting a <r + -polarized photon. The cylinder 
consists of silica, has a radius R = 250 nm and a refractive 
index n = 1.45. The y-axis is taken as quantization axis, 
(b) Schematic emission pattern of a circular dipole antenna 
close to a waveguide obtained from a finite-difference time- 
domain simulation. Shown is a heat-map plot of the electric 
field energy density. A part of the emitted light is coupled 
into the waveguide (green horizontal line, width of 100 nm) 
and shows a strongly directional behavior, as expected [9]. 
The emission into radiation modes is asymmetric as well. For 
simplicity, the simulation is performed in two dimensions. For 
illustrative purposes only, we assumed a refractive index of the 
waveguide of n = 3 > 1.45, which yields a larger asymmetry 
in the emission. 
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FIG. 2. (a) Partial decay rates 'y/2n into radiation and guided 
modes that propagate towards the +z or —z half space plotted 
as a function of the atom-surface distance (Ia . The horizon¬ 
tal black dashed line indicates the free-space limit Foo/47 t. 
(b) Directionality a of the total spontaneous emission. For 
parameters see text and Fig. [I] 


into a set of orthonormal guided and radiation modes of 
the nanofiber [251 . By solving the Schrodinger equation 
for the spontaneous emission problem [261 . we obtain a 
total emission rate that is the sum of partial decay rates 
into guided modes, and 'y^'] np into radiation modes. 
Here, the index p labels the polarization, / = ±1 indi¬ 
cates the propagation direction of the guided modes, k z 
is the projection of the wave vector of the radiation mode 
onto the fiber axis, and m € Z is the mode order. In gen¬ 
eral, the partial decay rates depend on va . We calculate 
the overall emission rates into the positive and negative 
half spaces as 


7± ?\r A ) = ^7/=±i, p (^), 

P 

^\r A ) = ±Y, 


±kl ° (R) ^ 

dk ^LL P (r A ), 


where the subscript indicates into which half space the 
emission is directed. 

The dependence of the partial decay rates on the atom- 
fiber distance (Ia is shown in Fig.[2](a). The guided modes 
are excited strongly asymmetrically, with emission to +z 
being about 10 times more likely than to —z. This is 
the result of the inherent link between local polariza¬ 
tion and propagation direction in strongly confined light 
fields [27J and matches the observations of a recent ex¬ 
periment jS], As the intensity of the guided nanofiber 


modes decreases with d A , so does the fraction of light 
that is coupled into them. The radiation modes show 
a complex behavior: The excitation of radiation modes 
that propagate into the +z or —z direction is, in gen¬ 
eral, also asymmetric. Moreover, the emission rates into 
these modes show characteristic Drexhage-type oscilla¬ 
tions [25] ■ Remarkably, depending on the radial distance 
of the emitter to the nanohber surface, either the radi¬ 
ation modes that propagate into the +z or into the — z 
direction are excited more strongly. The amplitude of the 
oscillatory behavior decreases with increasing emitter- 
surface distance. The radiative emission into the ±z 
half-space is symmetric for emitters far away from the 
nanohber, i.e. for {d A A 10 ). In this limit, the total 
emission rate of the atom approaches the free-space value 
of T 00 /(2 tt) = 5.234 MHz 

For an atom in the proximity of the nanohber, the com¬ 
bined emission rates into the ztz half space are in general 
not equal. The asymmetric emission should give rise to a 
force on the atom parallel to the fiber axis. We quantify 
the asymmetry by introducing a directionality parameter 
a, which is a sum over all partial decay rates weighted by 
the projection of the respective wave vector onto the hber 
axis. We denote the positive propagation constant of the 
guided modes at frequency u> 10 by /?io- The weights are 
then //?io/fcio for guided modes and k z /k\Q for radiation 
modes. Finally, we make a dimensionless by normalizing 
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to the total emission rate, 
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Figure [2jb) shows the directionality a as function of 
d A . An oscillatory, asymmetric emission is clearly visible. 
At the surface, the directionality reaches more than 20 %. 
The directionality is independent of the magnitude of 
the dipole moment. However, it strongly depends on the 
polarization of the emitted light: For the emission of a 
instead of a <r + -polarized photon, the directionality 
changes sign and, thus, can be controlled via the internal 
state of the atom. The emission is symmetric for the 
emission of 7r-polarized light. 

The described directional emission in conjunction with 
the conservation of total momentum in the system im¬ 
plies a lateral force on the atom that points towards the 
direction of stronger emission. This force can be viewed 
as a Casimir-Polder force on an atom at position r A . 
It is given by the ensemble-averaged Lorentz force [22j 
F = ■ E(r)}\ r of the quantized vacuum electric 

field E acting on the atomic dipole moment d. After 
solving the coupled atom-field dynamics and evaluating 
the averages over the atomic and field states, one finds 
for an atom initially in an energy eigenstate | n) a force 
F n (r A ) = F£ Ies (r A ) + F* es (r A ). The Green’s tensor can 
be decomposed into a vacuum part G (( ^ and a part G (1) 
associated to the scattering of the electric held from the 
nanofiber. Only the latter can give rise to a force. The 
non-resonant part of the force, F£ Tes (r A ), is proportional 
to the gradient of the symmetric part of G^ 1 '. The reso¬ 
nant part, F™ s (r A ), contains the gradient of the Hermi- 
tian part of G (1) [26l . 

When the atomic state |n) and the Hamiltonian are 
time-reversal symmetric, the dipole-matrix elements can 
be chosen real. In such cases, the gradient V can be 
replaced with a total derivative V rA , showing that no 
forces exist in a direction in which the system is transla- 
tionally invariant. We require an atom in an eigenstate 
that is not time-reversal symmetric, such that complex 
dipole-matrix elements can give rise to a lateral force. 

To see whether such forces can exist, let us con¬ 
sider an infinitely long cylinder. Its Green’s tensor is 
translationally invariant along the cylinder axis z , so 
that (d/dz)G l ' 1 \r 1 r',uj) = —(d/dz')G l - l \r,r',uj) holds. 
Combining this with the Onsager reciprocity relation 
G <1 - > (r,r',w) = G^^r', r, w), we find that the deriva¬ 
tive with respect to 2 : of the scattering part of the Green’s 
tensor is anti-symmetric. This immediately shows that 
for a ground-state atom with its purely off-resonant CP 
force a lateral force cannot exist. Such a lateral force is 
also forbidden by energy conservation: If it existed, one 
could use it to accelerate a ground-state atom along the 


fiber and, thus, gain kinetic energy while leaving the in¬ 
ternal energy of the atom and that of the environment 
unchanged. 

A lateral force can thus only arise due to the reso¬ 
nant component which is associated with the recoil of 
the atom when undergoing an optical transition between 
two states. These forces are fueled by the atom’s internal 
energy. Using the anti-symmetry mentioned above, the 
lateral force is given by 


F n,l( r A ) = 2i^o Y u 2 nk 

k<.n 

d 

X dnk ■ Im G (1) (r , r A , uj nk ) • d kn | r= ^ . (4) 

Being proportional to ImG, it is associated with the 
out-of-phase interaction of the electric-dipole oscillations 
with the reflected electric field. By contrast, the normal 
resonant CP force depends on the real part of the Green’s 
tensor. 

For the numerical evaluation of the force, we insert the 
scattering part of the Green’s tensor 


G ( 1 ) (r,r\o;) = ^ / dk z Y Y 

0 m =0 p,p'=TE,TM 


8n _ 

2 ^0m 


k 2 p 


Tpp'Clk z mp('r) < S)Q'—k z mp'('r ) ( 5 ) 


(<5,:y: Kronecker symbol) of a cylinder, which is given in 
terms of cylindrical vector wave functions a kzTnp and the 
respective reflection coefficients r pp > [30] with the disper¬ 
sion relation w 2 /c 2 = k 2 = fc 2 + fc 2 . Assuming a cesium 
atom in the excited state |?r) = |I) with a single down¬ 
ward transition to the ground state | k) = |0), and using 
the complex refractive index 1.45 + i2.05 x 10~ 7 of fused 
silica |3T] at the transition frequency wio, we can numer¬ 
ically evaluate the lateral force Q. We compare this to 
a calculation based on the solution of the Schrodinger 
equation where the expression for the lateral force based 
on the guided and radiation modes of the nanofiber is [261 


Fi, z ( r A) = ~h( f{ r a) 

V f,p 

/ &10 \ 

dkzk^l^rA )). ( 6 ) 

m,p '^ 10 ' 

Here, the recoil nature of the force becomes evident, as 
the summation is over partial decay rates into guided and 
radiation modes multiplied by the photon momentum of 
the respective mode. 

The result of the numerical evaluation of Eqs. @ 
and ([6]) is displayed in Fig. [3])a). One observes an os¬ 
cillating force with a distance dependence similar to the 
directionality parameter shown in Fig. [2jb) but which, as 
expected from momentum conservation, points into the 
direction opposite to the dominant emission. The lateral 
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FIG. 3. (a) Lateral Casimir-Polder force F\ iZ {dA,t = 0) on 
an excited cesium atom at distance as ck from the nanofiber 
surface, (b) Mean lateral velocity increase per emitted pho¬ 
ton. Force and velocity have been calculated using a method 
based on the solution of the Schrodinger equation, and using 
the Green’s tensor. Insets show magnifications for distances 
above 300 nm. We attribute the small deviation between the 
two methods to numerical artifacts. For parameters see text. 


force is an illustration of the subtle differences between 
Casimir-Polder force and potential [22]: It exists in spite 
of the ^-independence of the Casimir-Polder potential 
and hence cannot be derived from the latter. 

Recall that the lateral force is a result of the decay of 
the atom to its ground state and the associated emission 
of a photon. Thus, the time dependence of the ensemble- 
averaged force is given by Fi tZ (rA,t) = e~ rt Fi tZ (rA)- We 
thus obtain an average momentum kick per emitted pho¬ 
ton of 


A (p z ) = F^(r A )/T(r A ) = - ahk 10 , (7) 

where a is the directionality defined in Eq. ©■ Fig- 
ure §b) shows the velocity gained per photon, which 
can reach values of up to 0.8 mm/s in our situation. 

The agreement between both methods provides some 
understanding of the mechanism that underlies the lat¬ 
eral force. Note, however, that the Green’s tensor ap¬ 
proach is more general and also allows for the investiga¬ 
tion of situations where the fiber shows substantial losses 
that prevent the introduction of orthonormal modes. In¬ 
deed, an artificial increase of the imaginary part of the 
refractive index by five orders of magnitude to 2 x 10~ 2 
results in a 50 % larger lateral force, which is likely to be 
driven by the increased nonradiative decay that is auto¬ 
matically captured in the Green’s tensor formalism. 


To gain some further intuition into the lateral force, we 
consider the simpler case of an atom in front of a semi¬ 
infinite half space. In this case, the scattering Green 
tensor takes the form 

G (1) (r,r'.w) = — [ cl2fcl1 c ik«-(r-r')+ik x (z+F) 

8tt 2 J k x 

x ^ (8) 

p=TE,TM 


with polarization unit vectors e r7 ±, Fresnel reflection co¬ 
efficients r p and the dispersion relation w 2 /c 2 = k 2 = 

k 2 + (fcH) 2 . Substitution into Eq. Q then leads to a sim¬ 
ple expression for the lateral Casimir-Polder force in the 
retarded limit cLa 3> Any in this case, the wave-vector 
integral is dominated by the stationary-phase point at 
fell = 0 and we find for weak absorption: 


Fr z (dA) = 


Im((i 0 i ;Z dio ;a ;)cu 2 0 

8e 0 c 2 cP a 


sin(2eUfcio) r p (9) 


with r p — [n(wio) — l]/[n(tLi 10 ) + 1]. One sees that the 
lateral force is due to the interaction of the atom with its 
image dipole behind the surface of the half space. Also for 
this model case, we observe Drexhage-type oscillations. 

In summary, we have theoretically described a trans- 
lationally invariant lateral scattering force which arises 
from asymmetric spontaneous emission of a circular 
dipole emitter that is in close proximity to an optical 
nanofiber. In contrast to all lateral forces previously 
studied in Casimir and Casimir-Polder physics, the force 
does not rely on a corrugation of the material surface. 
Moreover, the magnitude and sign of the force depend 
on the polarization of the emitted light and thus can be 
controlled, e.g., by the quantum state of an atomic emit¬ 
ter. The described lateral force is generic in the sense 
that it prevails also in other geometries, e.g. for circular 
emitters above a plane surface. In contrast to the forces 
on Mie particles proposed in [3], our lateral force does 
not rely on a higher-order interaction between electric- 
and magnetic-induced dipoles. 

In optical force measurements on particles on a sub¬ 
strate [35], this effect will influence measurement out¬ 
comes as soon as scattering becomes relevant. Moreover, 
we expect the lateral force to enrich the dynamics of op¬ 
tically driven self-organization of atomic ensembles close 
to waveguides [231 !M!. In order to tune the force, the 
atom can be coupled to, e.g., a whispering-gallery mode 
optical resonator [351 3B], thereby changing the relative 
share of emission into guided and radiative modes. Lat¬ 
eral forces as described here will also influence laser cool¬ 
ing of atoms close to surfaces [37[ 38 j and nanophotonic 
structures PTTTj - PTT] . The study of lateral scattering forces 
might be extended to thermal Casimir Polder forces or to 
other geometries, e.g., an atom above a sphere where the 
lateral position is coupled to the atom-surface distance. 

We thank F. Ciccarello, F. Intravaia, Fam Le Kien, 
A. Rauschenbeutel, C. Sayrin, and J. Volz for helpful 






















5 


comments and discussions. Financial support by the 
DFG (grant no. SCHE 612/2-1) is gratefully acknowl¬ 
edged. S. Y. B. gratefully acknowledges support by the 
DFG (grant BU 1803/3-1) and the Freiburg Institute for 
Advanced Studies. 


Appendix A: Partial emission rates 


In this Supplement we present the calculation based 
on the solution of the Schrodinger equation for a cesium 
atom close to the nanofiber. The atom is initially in 
the hyperfine state |1) = | F' = 5 ,M' F = 5) of the 6P 3 / 2 
manifold. The only available decay channel is to the hy¬ 
perfine ground state |0) = |F = 4. M F =4) through the 
emission of a er + -polarized photon at frequency wio, and 
we can treat the atom as an effective two-level system. 
In the interaction picture, the atomic dipole operator is 
given by 


rotating-wave approximations is given by 
Hmt = -d • E(r A ) 

POO 

= -ifi V / dec G ufp (r A ) |1)(0| 

r J 0 

p£d/c 


f,p 


poo pcd/c 

— i h £/ dej / d k z G U }k zr np{f' a) 

J-uj/c 

X 11)(0 1 + H. c. 


(A4) 


We have introduced the coefficients G u f p (r) and 
Gukzmpi'r) which characterize the coupling of the atomic 
transition to a specific guided or radiation mode, respec¬ 
tively. These coefficients are given by 


G, 


Ulfp\ 


w/3' 
47T£o h 


dio ■ e w / p (r) e 


KfP z +p4>) 


G u k z m P {r) = 10 • e^fpir) e ^ z+m ^. 


(A5) 


d = doi |0)(1| e-^* + d w |1)(0| e iUlot , (Al) 

where doi = dj 0 is the corresponding dipole matrix el¬ 
ement. We follow [25] and decompose the electric field 
into contributions from guided and radiation modes, ne¬ 
glecting material absorption. For the guided modes we 
assume that the single-mode condition is satisfied for a 
finite bandwidth around wio- In cylindrical coordinates 
r = (r, <f>, z) , the positive-frequency part of the guided- 
mode field operator in the interaction picture can be writ¬ 
ten as 


4 +, W = i£ 

f,P 

X a w/p e w/p (r)e- i ^-^-^. (A2) 

Here, /? is the propagation constant of the guided mode, 
(J’ = d/?/dw, is the annihilation operator and e u f p 
is the profile function of the guided mode. The indices 
/ = ±1 and p = ±1 indicate the propagation direction 
and the handedness of the quasi-circular polarization, re¬ 
spectively [25] , Similarly, for the radiation modes 


dw 


huj/3' 

47re 0 


4 +) 


(r)=: 


m,p 1 


dw 


d k. 


' —UJ / C 


hui 

47Te 0 


x a wfc , TOp e u , fc , rop (r)e- i <" t -*«*- T "*>, (A3) 


For the combined atom-field state, we restrict the 
Hilbert space to a single excitation. In the following, 
the first element of the state vector indicates the state of 
the atom and the second element that of the field. Then, 
the atom-field state at any time can be written as 

pOO 

|^(t)) = Ci(t) |1,0) + V / dw c ufp {t) |0, l ufp ) 
f, P Jo 

POO PCD / C 

+ / dw / d k z Cuik^mp |0j lc dk z mp ) • (A6) 

m,p J 0 J-u/c 


We obtain expressions for the time-dependent coefficients 
c(t ) by inserting into the Schrodinger equation and fol¬ 
lowing a standard Wigner-Weisskopf treatment. First, 
we formally integrate the differential equations for the 
guided and radiation mode coefficients. Applying a 
Markov approximation, we replace c\(t') by Ci(t) in the 
resulting expressions and obtain 


Cufp(t) 0] (f) G 
C-l ok z mp{^) = Cl (t) G 


ufp(r A ) 

Jo 

Zk,mp(r A ) 

J 0 


(A7) 


These are inserted into the differential equation for c\ (t ), 
where we let the upper limits of the time integrals tend to 
infinity. Solving the integrals, we further neglect the term 
corresponding to a vacuum frequency shift and obtain a 
simple exponential decay, C\(t) = e~ Tt ^ 2 with a total, 
position-dependent spontaneous-emission rate r(r^) as 


where k z is the projection of the wave vector onto the 
fiber axis, m £ Z is the mode order, and p is the mode 
polarization. 

Let r A denote the position of the atom. The Hamil¬ 
tonian for the atom-field interaction in the dipole and 


^=£|G. 10/p (r A )| 2 

f,P 

rk 10 

+ £ / dkz \ G ^iok z mp{r A )\ 2 . (A8) 
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We see that each guided mode has an associated partial 
emission rate given by 


7 f P \ r A) = 27t \G Ul0 f p (r A )\“ 
w io/3io 


2 enh 


I d 


10 


ifp{ r A) | 


(A9) 


and similarly for the radiation modes, 


'Yk z mp( r - A ) — \G Ulo k z mp{ r A) | 
= 2cq/( ^ 10 ' 


(A10) 


To study the directional dependence of the sponta¬ 
neous emission, we note that the partial decay rates into 
each guided or radiation mode are proporti onal to the 
absolute square of the coupling coefficients (A5). The 


propagation direction along the fiber axis for the guided 
modes is encoded in the parameter /, such that / = +1 
corresponds to emission into the positive half space. Sim¬ 
ilarly, for the radiation inodes it is encoded into the pro¬ 
jection k z of the wave vector, and emission into the pos¬ 
itive half space is represented by modes with k z > 0 . 
Hence, the partial decay rates into the half spaces ±z for 
guided ( 7 ± G ' ) ) and radiation modes ( 7 )^) are given by 


(All) 


7± G) (n4) = 

v 

_ /*±fclO 

7± R) M = ± Y I dk * ^tlnp^A ) 

m,p ^ 0 

which is Eq. (2) in the main text. 


Appendix B: Lorentz Force 

The calculation of the Lorentz force amounts to com¬ 
puting the expectation value (F) = V(d • E{r))\ r=rA . 
We note that, in contrast to the interaction Hamilto¬ 


nian (A4), here the field operator is evaluated at a gen¬ 


eral position r. The position of the atom is only inserted 
after taking the gradient. 

We are only interested in the lateral part of the force, 
i.e., (F z ) = d ■ E(r))\ r=rA . To keep track of the 
individual contributions to the force, we label the part 
of the state ( |A 6 [ ) that describes excitations of the guided 
modes \4>g)i and similarly for the radiation modes \i/jr). 
The Lorentz force is then given by the sum of two terms, 
(F z ) = (F 2 (G) ) + (F 2 (R) ), where 


(FW) = 2Re(cl(t)-?-V liM (r) 


(Bl) 


and M £ {G,i?}. We have introduced the matrix ele¬ 
ment VqM = ( 1 , 0 | • E^\r) \iPm), where d^ is the 

negative-frequency part of the dipole operator, i.e., the 
second term in Eq. (Al). 


For the guided modes, we find 

POO 

Vi,c( r ) = -he ~ rt/2 Y / ^G ufp {r)Gl fp (r A ) 
f, P Jo 

e -i(u>-LU 10 )t _ 1 

X - 

U} - IV 10 

47re 0 

x y ] [dio • e UlBfp (r)][d 10 e uJio fp(fA)}* 

flP 

X e i[fPl0(z~ZA)+p((f>-(f>A)] 

(B2) 

where for the approximation we have made use of the fact 
that the mode functions, propagation constant, and their 
derivatives vary slowly in a frequency interval around 
Wio, for which we expect the main contribution to the 
integral in (B2). We hence replace them by their reso¬ 


nant values, such that the remaining integral is 


m = [ 

Jo 


a -i(LU-LU 10 )t _ ^ 


did U! - 

0 IV ~ LO IQ 


/ oo 

dA (wio + A) 

-OO 


0 -iA t 


- i 


(B3) 


= -17TW10. 


In the second line, we substituted A = iv — wio and ex¬ 
tended the lower limit of the integral from — wio to — oo. 
Inserting (|B3|) and (B2) into (Bl), we find that the lat¬ 


eral force due to spontaneous emission into guided modes 
is 

(FP) = £ / |d,o • e„„ /p (r,)|’ . 

zee — 


f,p 


_ (B4) 

The calculation for the contribution from the radiation 
modes is similar. The approximation that the main con¬ 
tribution stems from a narrow frequency interval around 
the resonance allows us to set the limits of the integral 
over k z to ±fcio and take it out of the frequency integral. 
Hence, 

V 1R (r)&- A ^ e Ft/ V [ 10 d 

X [^10 ’ ^ uJiok z mp ( r )][dl0 • &uj 10 k z mp{ r A)] (B5) 

such that 

., _ pk 10 

( F z R) ) = Y / dkz kz \ dw ' e u w k z mp{rA )\ 2 ■ 

Ze ° m,pj - k 10 

(B 6 ) 

Finally, we note that Eqs. (B4) and (B 6 ) can be rewrit¬ 


ten in terms of the partial decay rates of Eqs. (A9) 
and ( |A10P , 

(FP) = -he- Tt Yffho'y ( fi\r A ), 

f,p 

rk 10 ( B7 ) 

(FW) = - he~ rt Y / dk * k z ^l p (r A ). 


771,p 


1 —k 1 
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Appendix C: Resonant and non-resonant parts of 
the Lorentz force 


where S G = (G + G x )/2 is the symmetric part of a ten¬ 
sor. The resonant component is given by 


In general, the Lorentz force contains contribu¬ 
tions from resonant and non-resonant components, 
F n (r A ) = F£ Tes (r A )) + F* es (r A ). We give here their full 
analytical expressions based on the Green’s tensor G. In 
particular, the non-resonant component is [25] 




> t2 

Z Ukn 


x Vd nk • SG (1) (r,r A ,i£) • d fen | 


(Cl) 


F™(r A ) = 2„ 0 J2“ 2 nk 

k<n 

xVd nk -'HG il) (r 1 r A ,uj nk )-d kn \ r=rA , (C2) 


where FLG = (G + G* T )/2 is the Hermitian part. Here, 
to kn and d nk are the frequencies and matrix elements 
for electric-dipole transitions of the atom and G (1) is the 
scattering part of the classical Green’s tensor for the elec¬ 
tric field. 

The directionality parameter a given in the main 
manuscript can also be expressed in terms of the Green’s 
tensor as 




Im G xx + Im G zz 


-l 


(C3) 
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